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ADDENDUM 

RESEARCH DEPARTMENT - BRITISH BROADCASTING CORPORATION 
Mathematical Report No. F-025 (1965/44) 

GRAPHICAL SOLUTION OF CUBIC AND QUARTIC EQUATIONS 

1. INTRODUCTION 

The original Report F-025 gave not only graphical methods of obtaining 
reasonable approximations to the roots of cubic and quartic equations, but also dis- 
cussed iterative methods of improving the accuracy. The procedure suggested by 
Porter and Mack , however, appears to be superior, and this was not appreciated at the 
time of writing the original report. 

The procedure of Porter and Mack is here explained with numerical examples, 
in the hope that an engineer will thus be able confidently to solve a quartic equation 
to the degree of accuracy usually required in practice. Butler' 2 has given a clear 
account of more refined methods of achieving high accuracy. Ihe Addendum thus has no 
originality; its object is merely to make other people's work easier. 

The method can be extended to equations of higher degree, but attention here 
is confined to quartics, since the biggest 'hurdle' is usually to realise that a 
straightforward procedure exists at all for solving equations of degree above 2. 

The above-mentioned procedure in effect permits any equation of degree up to 
5 to be solved, since an equation of odd degree has at least one real root which can 
be found relatively easily (as indicated in Section 4 of the original report) and 
' factored out' . 

2. EXAMPLES OF THE METHOD OF PORTER AND MACK 

Consider the quartic equation 

X + Ci_X + Cq.X + c 3 x + Ca = (1) 



This can be factorized into 



(/ + ax + b) (x 2 + Ix + m) = (2) 



x A + (a + l)x 3 + (6 + al + m)x 2 + (61 + am)x + bm = (3) 



2. 



Equating coefficients of Equation (3) with those of the original quartic (1) we have 
the following relations 

a + I = Cl (4) 

b + al + m = c 2 (5) 

6Z + am = c 3 (6) 

6m = c 4 (7) 

The method of Porter and Mack consists in choosing a likely value b for b, 
and then solving Equations (4), (6) and (7) for the corresponding values of a, I, m 
giving 

m r = cjb r (8) 

a r = C -p-^l (9) 
o - m. 



Z = Cl - a (10) 

r *■ r 

where a , I , m are the corresponding values of a, Z, m when 6=6 respectively. 

Equation (2) is then written as 

g =6 + a I. + m - c 2 (11) 

We seek to choose b so that g is zero. In other words, we choose 6 , calculate a , 
Z and m as if this value were correct, and then the value of g tells us how close 
we are to the correct value of b. 

Consider first the case when c* is positive. Test to see whether c^ - 
(cg/ci) . This is the only 'special' case, as Butler has pointed out; we then have 



while a and Z are found from 



so that 



6 = m = Vca. (12) 

a + I = ci (13) 

6 + al + m = c 2 (14) 

aZ = c 2 - 2/ca (15) 
a - Z = ( Cl 2 - 4aZ)* 

= {ci 2 - 4c 2 + 8/c 4 }* (16) 

a = J4(ci + [ Cl 2 - 4c 2 + 8/c 4 ]*} 

Z = ci - a (17) 



Example 1 



x A + 5/ + 8% 2 + 5x + 1 = (18) 



2 

Ca is positive and equal to (c 3 /ci) 

b - m = Vc 4 = 1 (19) 

a = V 2 {5 + (25 - 32 + 8)*} = 3 (20) 

I = 5 - 3 = 2 (21) 

2 2 \ 

giving (x + 3x + 1) (% + 2% + 1) I 

> (22) 

= x 4 + 5x S + 8* 2 + 5* + 1 ) 

In general, however, c* 7^ c 3 /ci and there is a positive value of 6 between 

and /c 4 . 

Take as first trial value of 6, ! ^ /c 4 . Calculate m, a, i and g from (8), 
(9) and (10). 

If g is positive, take next value of 6 as % /c 4 . 

If g is negative, take next value of b as % vc 4 . 

Calculate m, a, I and g again from (8), (9) and (10). If both values of g are 
positive, take next value of b as O'Wc^.. 

If both values of g are negative, take next b as O'lvc^. 

If the 2 values of g are of opposite sign, then the next b to try is calculated by 
linear interpolation, i.e. 

b 3 = blg * ' b2gl (23) 

g2 - gl 

where bi, b 2 and b 3 are the successive b's, g t and g 2 being the g's associated with 
bi, b 2 . 

Then calculate m, a, I and g from (8), (9) and (10) with 6 = b 3 . 

If the g's are still remaining either all positive or all negative, the b's to be 
tried are respectively Q'99/ca, 0'999v / c 4 etc. and O'Ol/ct, O'OOk/c*. etc. Otherwise, 
linear interpolation is used to find the successive b's. The process in continued 
until g is small enough for the accuracy required. 



Example 2 



% 4 + 2/ + 10/ + Sx + 1 = (23) 



C4. is positive, and not equal to c 3 /ci so we try first 

6 = x A^c A = 0'5 
The working is best set out as in Table 1. 

TABLE 1 

c± = 2, c 2 = 10, c 3 = 5, C4 = 1 



6 

r 


m 
r 


b - m 
r r 


cib r - c s 


a 

r 


l r 


S r 


0-5 


2 


-1-5 


-4 


2-67 


-0-67 


-9-29 


0-25 


4 


-3-75 


-4-5 


1-2 


0-8 


-4-79 


0-1 


10 


-9-9 


-4-8 


0-485 


1-515 


0-835 


0-12 


8-333 


-8-213 


-4-76 


0-5796 


1-420 


-0-724 


0-111 


9-009009 


-8-898009 


-4-778 


0-536974 


1-463026 


-0-094384 


0-1097 


9-115770 


-9*006070 


-4-7806 


0-530820 


1-469180 


0-005340 


0-109770 


9-109957 


-9-000187 


-4-78046 


0-531151 


1-468849 


-0-000092 



This gives for the 2 quadratic factors of the left hand side of (23) 

(/ + 0-531151* + 0-109770) (* 2 + 1-468849* + 9-109957) 
= * 4 + 2-000000* 3 + 9-999908* 2 + 4-999998* + 1-000000 



(24) 



It should be noted that there is no need to work to many significant figures until g 
becomes small. 



Example 3 



usual) . 



(Illustrating the case when the g's remain positive for longer than is 



x A + 5x S + 8-01x 2 + 5-03* + 1-01 = 



(25) 



c A is positive, and nearly but not exactly equal to c 3 /ci . 
Try first, b = Vi^c* - 0*5. The working is as in Table 2. 



TABLE 2 
Ci =5, c 2 = 8-01, c 3 = 5-03, c 4 = 1*01 



b 

r 


m 

r 


6 - m 

r r 


cib r - c 3 


a 

r 


I 

r 


g r 


0-5 


2 


-1-5 


-2-5 


1-7 


3-3 


0-1 


0-75 


1-65 


-0-6 


-1-28 


2-13 


2-87 


0-203 


0-91 


1-1099 


-0-1999 


-0-48 


2-4012 


2-5988 


0-2501 


0-995 


1*015 


-0-020 


^0-055 


2-75 


2-25 


0-1875 


1-0040 


1-0060 


-0*0020 


-o-oioo 


5-00 





-6-00 



The last entry in Table 2 indicates that g is varying violently in this region because 
b m m . Linear interpolation may not therefore be satisfactory as a means of obtain- 
ing the next b , but the correct b must be between 0*995 and 1*004. As an alternative 
to linear interpolation, therefore, the next 6 can always be taken as the average of 
any two 6 associated with g of opposite sign, and in this way the interval within 
which b must lie can be progressively reduced until eventually linear interpolation is 
adequate for the degree of accuracy required. In this particular case the required 
quadratic factors are exactly (x + 3x + 1) (x + 2x + 1*01). 

Next consider the case when c 4 is negative. We first determine 

H = ( Cl 2 /4) + (c 3 2 /4c 4 ) - c 2 (26) 

If H>0, there is a value of b between and vj c^. | . Take as first value of b, Vi v | c^. | . 

Calculate m, a, I and g from (8), (9) and (10) as usual. 

If g is positive, take next value of b as % v\ct r \. 

If g is negative, take next value of b as % v\c A \. 

Recalculate m, a, I and g. 

If both values of g are positive, take next value of 6 as * lv | c 4. | . 

If both values are negative, take next 6 as * 9v | C4. | . 

If the two values are of opposite sign, the next b is obtained by linear interpolation. 

Recalculate m, a, I and g. 

If the g's are still remaining all positive, the b ' s to be tried are 0-0l/|c 4 |, 
0-00l/|c4.| etc. Otherwise, linear interpolation can be used, or the alternative 
averaging procedure mentioned in Example 3. 



Example 4 



x 4 + 9x 3 + 2x 2 + 5x - 2 = (27) 



So we try, 



c 4 is negative. H = 15*125>0 

b = ViV\cA = 0-7 

ci = 9, c 2 = 2, c 3 =5, c 4 = -2 



T/lBLtf 3 



c\b -c c 



Z 



8*63 -0*97 

7*459 8-667 

8-518 0-143 



0-7 -2-86 3-56 1-3 0*37 

1-06 -1-887 2-947 4'54 1-541 

0-74 -2-703 3-443 . 1'66 0'482 

0-735 -2-7211 3-4561 1-615 0-46729 8*5327 0-00115 

0-73496 -2-72124 3-45620 1*61464 0*467172 8*53283 0-000019 

which gives for the two quadratic factors required 

U 2 + 0-467172* + 0-73496) (/ + 8 -53283% - 2-72124) 
= / + 9* 000002a; 3 + 2*000019% 2 + 5-000002* - 2-000003 



(28) 

If H is negative, there is a value of b between V|c 4 | and 0. Take as first value of 

b, -%Ac*\. 

Calculate m, a, I and g using (8), (9) and (10) as usual. 

If g is positive, take next value of b as -% v|c 4 |. 

If g is negative, take next value of 6 as -M v|c4.|. 

Recalculate m, a, I and g. 

If both values of g are positive, take next value of b as -0" 9k | c 4 | . 

If both values of g are negative, take next b as -0'lv |c 4 |. 

If the 2 values are of opposite sign, the next b is obtained by linear interpolation. 

Recalculate m, a, I and g. 

If the g's are still remaining all negative, the 6's to be tried are -0-0l/|c 4 |, 
-0*001v|c 4 | etc. respectively. Otherwise linear interpolation is used, or the 
averaging procedure mentioned in Example 3. 
Example 5 



/ - 5x 3 + 7/ + 9x - 4 = 



(29) 



So we try 



25 81 

C4 is negative. H 3 7<0 

4 16 



b = -W\c A \ = -1 

ci == -5, c 2 =7, c 3 = 9, c 4 



T/4ELJ? 4 



6 

r 


m 
r 


b - m 

r r 


c 1 b r -c 3 


a 

r 


I 

r 


g r 


-1 


4 


-5 


-4 


0-8 


-5-8 


-8-64 


-0-5 


8 


-8-5 


-6-5 


0-76 


-5-76 


-3-88 


-0-2 


20 


-20-2 


-8 


0-396 


-5-396 


10-6632 


-0-42 


9-524 


-9-944 


-6-9 


0-694 


-5-694 


-1-8476 


-0-388 


10-30927 


-10-69727 


-7-06 


0-659981 


-5-659981 


-0-81421 


-0-3628 


11-025358 


-11-388158 


-7-186 


0-631006 


-5-631006 


0-109359 


-0-36578 


10-935535 


-11-301315 


-7*1711 


0-634536 


-5-634536 


-0-005561 


-0-365636 


10-939841 


-11-305477 


-7-17182 


0-634366 


-5-634366 


-0-000045 



This gives for the 2 quadratic factors 

(% 2 + 0-634366* - 0-365636) (/ - 5-634366* + 10-939841) 

= / - 5- 000000* 3 + 6- 999955/ + 8-999990x - 3-9999997 



(30) 



If H =0, we can take b = -1/ | c 4 1 , m '-= +(/ 1 c 4 1 , a and I are then found from Equations (4) 
and (5) which give 



Example 6 



a = V 2 { Cl + ( Cl 2 - 4c 2 )*} 
/, - c x - a 



/ + 6x 3 + 7/ + 4* - 2 = 



c 4 <0 H = 

6 = V2 = -1-414214 

m ~= +1-414214 

a = 4-414214 ■ 

1=6- 4-414214 = 1-585786 



(31) 



(32) 



giving 



(/ + 4-414214* - 1-414214) (/ + 1*585786* + 1 '414214) 
= x 4 + 6-000000/ + 6-999999/ + 4-000002* - 2-000001 (33) 
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November 1965 Mathematical Report No. F-025 

(1965/44) 

GRAPHICAL SOLUTION OF CUBIC AND QUARTIC EQUATIONS 



SUMMARY 

Given any cubic or quartic equation with real coefficients, it can be 
transformed by a change of variable into a standard form, and thus s'olved by means 
of a double- sided scale (Fig. 1) in the case of a cubic, or by means of the nomogram 
of Fig. 3 (or Fig. 4) in the case of a quartic. Alternative procedures for certain 
cases where high accuracy with this nomogram is difficult to obtain are also discussed. 
Thus the information here given enables cubic and quartic equations to be solved as 
easily (to an accuracy of say two significant figures) as quadratic equations are 
solved without this information. Once a good approximation to the roots is known, 
accuracy can be improved by iterative processes if necessary; this is briefly dis- 
cussed. The normal routine required for solving a cubic or quartic to graphical 
accuracy is specified in the Conclusions (Section 5). 



1. INTRODUCTION 

Various well-known methods are available for solving algebraic equations 
with real coefficients, including special ones which are available only for cubic or 
quartic equations. Most of these involve either iterative procedures (which can be 
tedious if no approximate prior information about the roots is available) or some 
transformation of the given equation into a special form, containing parameters which 
have to be determined by means of an auxiliary calculation. It is reasonably easy to 
obtain approximate values of real roots by drawing a graph of the left hand side of the 
equation, but in the case of an equation of even degree all of whose roots are complex, 
there is no obvious method of locating complex roots approximately except the now well- 
known 'root-locus' method which again is straight -forward in principle, but usually 
requires a considerable amount of calculation before an adequate approximation to the 
roots being sought becomes available. 



°-6 CP 
0025 003 0-04 0-05 06 C-1 0-2 0-3 *l 1 2 3 4 5 55 6 63 

| 1 I I ,1 I L-jJ • 1— J I | I I pi I I [ 1—1-1 1 1 I I 1 1 1 1 1 1 1 I I I I I I l_ I | I I Lpl J !_,_! I I ■ ■ ■ ■ | ■ ■ ■ ■ I ■ ■ I 
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Fig. 1 - Double-sided scale for finding the Positive Real Root of Equation (5) 
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Fig. 2 - Nomogram for finding the number f3 in expression (18) associated 
with the factors of the left hand side of equation (16) 



The general method of attack here considered is to use a two-sided scale 
(Fig. 1) to obtain the real root of a cubic, and a nomogram (Fig. 2, 3 or 4) to obtain 
the number-term of a real quadratic factor of a quartic. This is far from being the 
first time that nomograms have been used for solving algebraic equations, since 
d'Ocagne • '^ gave techniques as long ago as 1902 for finding the real roots of equa- 
tions of degrees up to 7, and Mack and Tomlin constructed a nomogram for the solution 
of quartics in 1945. 



Any algebraic equation in x can be 'normalized' by changes of variable so 
that there is no term in x n ' % if the equation is of degree n, and so that the term 
independent of x is ± 1. This normalization must be carried out before the techniques 
to be described are used. If the equation so normalized is cubic, the relation 
between the positive real root and the normalized coefficient of x is one that can be 
determined and expressed graphically as in Fig. 1 once for all as is explained in 
Section 2. If the equation so normalized is quartic", the relation between the 
coefficient of x , the coefficient of x and the number term in the quadratic factor 
being sought can be expressed in the form of a nomogram (Figs. 2, 3 or 4) as is 
explained in Section 3. Cases where there is at first sight some difficulty in using 
the nomogram satisfactorily are also discussed in Section 3. When a good approxima- 
tion to a root of an equation is known, accuracy can be improved by means of iterative 




Fig. 3 - Nomogram analogous to Fig. 2 for the case when b Q in equation (Ik) is negative 

procedures; this is briefly discussed in Section 4 and more fully in References 4, 5 
and 6. Section 5 gives conclusions, in the form of a specification of the routine to 
be followed for solving a cubic or quartic. 



2. CUBIC EQUATIONS 

Suppose that the given equation is 

then by writing 



3 2 

x + a^x + aix + a = 



y ~ x + —a 2 



Equation (1) can be reduced to 



y + b lV + b = 



where 



. 1 2 1 2 3 

bi = a t - -a 2 ; b = a o --a.!a 2 + r-a 2 



(1) 



(2) 



(3) 



(4) 



C,.( 
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Fig. 4 - Nomogram identical with Fig. 2 except for the suppression 
of the negative-/3 branch 

l/s 
and finally y can be changed to -6 ' z so that the equation reduces to the normalized 

form 



:i z - 1 = 



(5) 



where 



ci = -6i/6_ 







(6) 



If 6 W 0, x fa - a 2 /3 is the required approximation to a real root of (1). The 
remaining factor of the left hand side of (1) can then be taken as 

2 2 3a„ 
* + -a 2 * + —2 
3 a 2 



Equation (5) necessarily has a positive root say a; the left-hand side is therefore 
expressible in the form 



2 1 

(z - a) (z + az + — ) 
a 



(7) 



and 



2 1 

ci = a - - (8) 

a 



Now it is easy to construct as in Fig. 1 a two-sided scale in which there is a linear 
Ci-scale on one side, and a non-linear scale of the corresponding values of a on the 
other side. From this double scale, the value of a corresponding to a given value of 
C\ can be read, and the factors (7) of the left-hand side of the normalized equation 
(5) deduced immediately. It only remains to express factors in terms of the original 
variable x; to do this we must substitute 

z = - U+ia^/b// 3 (9) 

into (7). The factors of the left hand side of the original equation (1) thus becomes 
respectively 

1 1/3 
x + — a 2 + 6 ' a = x + p, say (10) 



and 



2 ' ' 2 ■ l/ ^+ga 2 "--^a6 1 / a + (6 o a / a /a)] 

2 
~ x + qx + r, say (11) 



+ x(— a 2 - aJb 
3 



Examples: 

In this and later sets of examples, cases with exact roots are chosen and 
the exact roots are given, so that the accuracy obtainable by the graphical technique 
is easily assessed. 

(i) x 3 + 9x 2 + 20x + 12 = 

or (x +1) (x + 2) (x + 6) = 

(ii) x S + 7x 2 + Ux + 8 = 

or (x + 1) (x + 2) (x + 4) = 

(iii) x 3 + 6* 2 + 11* + 6 = 

or (x +1) (x + 2) (x + 3) = 

(iv) x 3 + 3x 2 + 4'25x +3=0 

or (x + 1*5) (x 2 + l'5x + 2) = 



Example No. 


a? 


cii 


a 



61 


fc o 


Cl 


a 


P 


9 


r 


(i) 


9 


20 


12 


-7 


6 


2-12 


1-56 


5-84 


3-16 


2-61 


(ii) 


7 


14 


8 


-2-333 


0'740 


2-86 


1*9 


4*05 


2-95 


1-865 


(iii) 


6 


11 


6 


-1 





- 


- 


9 


4 


3 


(iv) 


3 


4-25 


3 


1-25 


0*75 


-1*51 


0*6 


1-545 


1-455 


1-83 



3. QUARTIC EQUATIONS 

Suppose that the given quartic equation is 

4- 3 2 

x + a 3 x + a 2 x + a t x + a = (12) 

then by writing 

1 
y = x + T a 3 (13) 

4 

it can be reduced to 

y A + b a y* + b iy + b Q = (14) 



2 

62 = a 2 ~ 0'375a s 

61 = at - 0"5a 2 a s + 0*125a 3 ^ (15) 

6 o = a o - 0'25a!a 3 + 0'0625a 2 a 3 2 - 0-01171875a 3 4 

Consider first the case in which b is positive, and put y - (b ' )z, and (14) 
reduces to • 

2 4 + C 2 Z 2 + CiZ + 1 = (16) 

where 

c 2 - 6 2 /b *; d - d/6 s/4 (17) 

We seek to factorize (16) - the factors will be of the form 

I 
(z 2 + az + p) (z 2 - az + -) (18) 

3 2 

so that the z term and the term independent of z are correct; in order that the z 

and z terms of (16) may be also correct, we must have 

„ 1 2 
£ + - - a = c 2 (19) 

- a(/3 - -) - Cl (20) 

Eliminating a from (19) and (20) we have 

2 

c 2 (21) 






Now Equation (21) is of the form 

hh + /2g3 + h a - (22) 



where 

2 



/i = Ci ; / 2 = c 2 ; / 3 = 1/0-C1//3)]- 
g 3 = i; h a =- [/3 + (1/0)] 



(23) 



so that /i depends only on Ci, / 2 depends only on c 2 , and f 3 , g 3 , h 3 depend only on /3. 
Equation (22) is a standard form, which means that the relation can be expressed by 
means of a nomogram (see reference 7) - it is the condition that the points 

P (-d o , /A/i), Q(d o , /Xj/ 2 ) and Rtf 3 , Y a ) 

shall be in a straight line, where d Q , /^ and fi 2 are scale factors at our choice, and 

j s = d . y g = (24) 

Migs + M2/3 Migs + M2/3 

As ci varies, the point P just mentioned moves along a line parallel to 
the 7-axis, and the lowest point of the scale occurs for ci = 0, since fj_ ~ c\ . On 
the other hand, as c 2 varies, the point Q also moves parallel to the F-axis, but the 
point corresponding to c 2 = should be near the middle of this scale since / 2 is 
proportional to c 2 and may therefore have either sign. There are many possible 
choices of d , \±\ and ji?; it is convenient to have d (half the distance between the 
point representing c\ - and the point representing c 2 = 0) a round number, and to 
decide in advance the approximate overall size of the nomogram. Fig. 2 was therefore 
constructed originally to have an overall size of 240 mm by 180 mm, and with the zero 
of the c 2 - seale 100 mm above the lower edge of the paper, so that 



d = !4/(240 2 + 100 2 ) = 130 mm; 



Mi = M2 - 2 mm per unit of corresponding / 



(25) 



The J-axis is then the oblique line joining ci - (bottom left-hand corner) to c 2 = 
(about half way up the right hand side of the diagram) and the Y - axis is vertical. 
The origin is midway between the vertical ci - scale on the extreme left and the 
vertical (and linear) c 2 - scale on the extreme right. The /3-scale, defined by 
equations (24), is curved, and has two branches, according to whether /3 is positive or 
negative. The0-scale has been calibrated with a value of /3 (numerically greater than 
unity). Once a value of has been found, Equation (20) gives the associated value 
of a and completes the factorization (18) to graphical accuracy. 

If ci is small, it may not be easy to determine accurately the position of 
the intersection with the /3-curve of the line joining the points (~d , /J-ifi) calibrated 
by the value of c± and the point (d Q , j-Lzf?) calibrated by the value of c 2 . When this 
happens, or if the index line meets the /3-curve beyond the last calibration /3 = +10, a 
starting approximation can be obtained by neglecting c 1( in which case (16) becomes 
quadratic in z so that 

1 

+ ~ «c 2 (26) 



Having found j3 from (26) a first approximation to a can be obtained from (20) and an 
improved value of /S is then obtainable from (19) and an improved (and probably ade- 
quate) value of a from (20). 

Equation (26) is unsatisfactory if |c 2 | <2. Instead, replace j3 initially 
by 1 in. (19) if 0<c 2 < 2 or by -1 if 0>c 2 > -2; this gives a first approximation to a, 
and (20) then gives an improved value of /3 which in turn can be substituted into (19) 
to give an improved value of a and so on. Alternatively, y can be replaced by c 2 as 
an approximate root of (36) below. 

If ci = and c 2 is greater than 2, the quartic (16) has four purely imagi- 
nary roots; we can therefore expect (16) to have four complex roots when ci is suffi- 
ciently small and c 2 >2. In Fig. 2 the index line (joining the point marked c 2 to the 
point marked c 2 ) in such a case intersects the positive-/3 branch of the /3-curve in one 
point and does not meet the negative-/3 branch at all. The factor (z + az + jS) thus 
derived will be associated with a complex root-pair of (16); the remaining factor 
z - az + (1//3) may be associated with two real roots of the same sign, or may be 
associated with a complex root-pair. If however c 2 is sufficiently negative, the 
possibility arises that the index line may meet the negative-/3 branch of the /3-curve 
in two places as well as the positive-/? branch in one. This is the case when (16) has 
four real roots, two positive and two negative, which can be combined in three dif- 
ferent pairs as real quadratic factors of the left-hand side of (16). 

As before, it is necessary to express the factors found in terms of the 
original variable x; this can be done by means of the substitution 

z = (x +la 3 )/b 1 ^ (27) 

4 ° 

and the real quadratic factors of the left hand side of (12) thus become 
x 2 + ( l Aa s + db^hx + (~a 3 2 +-a 3 a6„ 1 / 4 + #> *) 



3 ' — u 3 u.u o ■ ■ yu Q ) (28) 



2 



~ x + pix + q lt say 
and 



x + (^a 3 - ai// 4 )* + (-i a 3 2 --agab// 4 + [&„*/#) (29) 

16 4 



x + p 2 x + q 2 , say 



Examples [6 Q >o] 



(v) x A + 4* 3 + 9x 2 + 8x + 5 = 

or (x 2 + x + 1) U 2 + 3x + 5) = 
(vi) x 4 + 3-2/ + 4-71* 2 + 3-04* + 0-8^0 

or (/ + l'l* + 0-4) (/ + 2'lx + 2) = 
(vii) x 4 + 8x 3 + 22-75/ + 29-5% + 17-5 = 



or (x + 2-5) (x + 3-5) (/ + 2x + 2) = 
(viii) / + 10* 3 + 35* 2 + 50% + 24 = 

or (* + 1) (x + 2) (* + 3) (* + 4) = 

Example No. 6 2 bi b c 2 Ci /? a p t g^ p 2 g 2 

(v) 3-2 3 1-732 -0-877 1-75 0-744 2-98 5-01 1-02 1-01 

(vi) 0-87 -0-4 0-1536 2-220 -1-630 2-1 1-00 0-97 0-33 2-23 1-96 

(vii) -1-25 2-5 1-5 -1-021 1-845 1-54 -2-072 1-71 1-30 6-29 9-38 

(viii) -2-5 0-5625 -3.1/3 0-3 5 4 5 6 

If 6 is negative, the relation between y and z is y = (-b ) ' z, and (14) 



reduces to 



z 4 + az 2 + ciz - 1 = (30) 



,3/4 



where 

c 2 = 6 2 /(-6 ) s ; Cl = bj{-b )°l* (31) 

and the quadratic factors are of the form 

(z 2 + az + f3) (z 2 - az - -) (32) 



We then have 

1 

- a(/3+ 1) = Cl (34) 

so that the relation on which the nomogram for this case is based (obtained by elimi- 
nating a between (33) and (34)) is 

08--)- Cl2 



/3 a = c 2 (33) 



i,,.; i* c * (35) 



M 



This is as before of the form (22) but 



/i = Cl 2 ; / 2 = c 2 ; / s = 1/03 + (1//3)] 2 (36) 

g 9 = 1 /i 3 = U/8) - /3 . 

and with the definitions in Equations (36) instead of Equations (23), Equations (24) 
and (25) still apply to define the /3-curve for the new nomogram (Fig. 3); the Ci 
and c 2 scales are unchanged. In Fig. 3, as in Fig. 2, /3 is specified as numerically 
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greater than 1; j3 = ±1 on the A"-axis, and (3 is positive above this axis and negative 
below it. If there are three intersections of the index line with the /3-curve, the 
quartic (30) has four real roots. If the index line meets the /3-curve beyond f3 = ±10 
the approximate value of /3 can be obtained from (33) with a - 0. 

Having determined j3 from Fig. 3 and deduced a from (34), we again have to 
express the factors found in terms of the original variable x; this requires the 
substitution 

z = (x + - 1 a s )/(-6 ) l/4 (37) 

*** 4 

and the real quadratic factor of the left hand side of (12) thus become 



x 2 + (&i a + at-bj 1 / 4 )* + |"ia s 2 + 'a 3 a('b )^ A + A -*>„)*] 



(38) 



2 

x + Pi_x + Qi say 



x 2 + (%a a 



il-bXh x + R-as 2 - K.ai-bJ 1 ^ - (-b o h//3] (39) 



2 

~ x + P 2 x + Qz, say 



Examples b <o 



(ix) x 4 + 3-2x B + 4'41% 2 + 2-41% +0-2=0 

or (x 2 + 1-lx + 0-1) (/ + 2-lx + 2) = 

(x) x 4 + 8x 3 + 19x 2 + 22x + 10 = 

or (/ + 6x + 5) (/ + 2x + 2) = 

(xi) x* - 25x 2 - 60x - 36 = 



or (x +1) (x + 2) (x + 3) (x - 6) = 



(xii) x 4 + 8x s + 27x 2 + 70x + 50 = 



or (x + 1) (x + 5) (x 2 + 2x + 10) = 



Example 6 2 6i 6 Q c 2 ci /3 a P a Q x P? () 2 
No. 

(ix) 0-57 -0*55 -0*1344 1*555 -2-478 2*0 0*991 2*2 1*85 1*0 -0*024 

(x) -5 10 -6 -2*041 2*608 -2*0 1*043 5*63 2*37 2*37 1*96 

(xi) -25 -60 -36 -4-167 -4*082 1-0 2-041 5 6-5-6 

(xii) 3 26 -30 0-054 2-178 1*4 -1-03 1-60 6-8 6-40 4-91 
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4. IMPROVING THE ACCURACY OF A GRAPHICAL FIRST APPROXIMATION 
RY ITERATIVE PROCEDURES 

The main difficulty in seeking complex roots of algebraic equations is to 
obtain a good starting approximation, and this we have succeeded in doing for cubics 
and quartics by the procedure of Sections 2 and 3. If we want greater accuracy than 
is available by graphical means, two alternative iterative methods of improving accu- 
racy are briefly discussed below. 

For a cubic equation, it is naturally the real root which should be found to 
a sufficient degree of accuracy in this way, and this root should then be factorized 
out, leaving only a residual quadratic to be solved in the usual way. If the real 
root has smaller modulus then either of the remaining roots, the processes to be 
described will be quite straightforward. If it has larger modulus then the remaining 
roots, it may be necessary to change z into 1/t in Equation (5), and actually solve 
the t-equation instead of the z-equation. 

Suppose that zftja is the approximation already obtained graphically and that 
the accurate value of z is a + S. Then we can write (5) in the form 



(z - a - S) |z 2 + (a + S) * 



(40) 



a + S_ 
and choose S so that the coefficient of z in the product is correct, that is, 

Cl = (a + §) 2 -Sa ! - (1/a) + S [2a + (1/a 2 )] (41) 

a + b 

If a is a good approximation, ci^sa - (1/a) and thus 8 is small; the process can be 
repeated with (a + S) instead of a if necessary. This procedure is effectively an 
adaptation of Rairstow's method. An alternative process originally due to Lin is 

even simpler in principle; we divide the left hand side of (5) by z(z-a); the remain- 
der is 

(a 2 - ci) z + 1 (42) 

and the next approximation to z is l/(ci - a ) = a say. If the process is now 
repeated and a further approximation a" thus obtained, then a very close approximation 
is given by Aitken's formula 

#2 /i .2 

a - a a A j 



= a + — = A (43) 



2a' - a - a" 2Aj. - A, 

where 

Aj = a' - a A2 = a" - a (44) 

If necessary, the whole process can be repeated with A as the starting approximation 
instead of a. 

When we are similarly trying to improve the approximation to a quadratic 
factor z + az + f3 of (16) or (30) the Lin procedure can be applied, not to (16) or 
(30) themselves, but to (21) or (35). In (21) put 



12 

+ (1//S) - 7 (45) 



and (21) becomes 



y(y 2 _ 4) _ Cs (y 2 _ 4 ) _ Cl 2 = Q 



7 3 - c 2 7 2 - 47 + (4c 2 - d 2 ) - (46) 



We have an approximation J3 to /3 from Fig. 2 or 4, and from this we can derive y = 
fi Q + (1//3 ) which is approximately a root of (46). We can therefore divide the left 
hand side (46) by 7(7 - y o ) and find a better approximation analogous to A in (43) and 
thus obtain 7 as accurately as we need. The corresponding value of /3 is obtained from 
the quadratic equation (45). In (35) put 



and (35) becomes 



fi- (1//3) =7 (47) 



y(y 2 + 4) - c 2 (y 2 + 4) - Cl 2 



7 3 - c 2 7 2 + 47 - (4c 2 + ci 2 ) = (48) 



and if /3 is an approximation to /3 obtained from Fig. 3, 7 ~ fl ~ (l//^.) is an 
approximation to a root of (48), so the left hand side can be divided by 7(7 - 7 ) 
and a better approximation analogous to A in (43) can be found as before. Equation 
(47) gives the corresponding value of /3. 



5. CONCLUSIONS 

The following are the normal routine steps required to solve any cubic or 
quartic equation 'to graphical accuracy'. 

A. Cubic Equations 

(i) Convert the given equation (1) into the normalized form (5), using equations 
(4) and (6). 

(ii) From Fig. 1 given a. from equations (6) and (4), read off a, a positive real 
root of (5). 

(iii) Substitute into (10), that is find p in order to obtain the linear factor of 
the left hand side of (1), and into (11), that is, finding q and r in order 
to obtain the complementary quadratic factor, and check against (1) by multi- 
plying out. The roots are then x = -p and x - -q + (q - 4r) 2 . 
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B. Quartic Equations 

(i) Convert the given equation (12) to the normalized form (14) by using equa- 
tions (15) and again to the form (16) using equations (17) if 6 >0 or (30) 
using equations (31) if b <0. 

(ii) Use the nomogram (Fig. 4 if b o >0, Fig. 3 if b <0) to determine a value of /3, 
the number term of the quadratic factor (18) or (32) being sought. (In 
Fig. 4 the negative-/^ part of Fig. 2 has been suppressed; this merely tells 
us whether the quartic has four real roots or not, and by factorizing the 
quadratics (18) obtained with positive /3, this information is equally well 
available. ) 

(iii) Use Equation (20) or (34) to determine a, the coefficient of z in the quad- 
ratic factor (18) or (32) being sought. 

(iv) Re-convert to the original variable x by substitution into (28) and (29) 
(b >0) or into (38) and (39) (b <0) and check against (12) by multiplying 
out the converted factors. 

If it is difficult to determine the intersection of the index-line with the 
/3-curve accurately enough (when c± and c 2 are small and 6 >0) take j3 ~ 1 as the first 
approximation if 0^02^2, /3 = -1 as the first approximation if 0^c 2 > -2 and use (26) to 
obtain the first approximation to fi otherwise. Then use (20) to obtain a, and (19) 
and (20) alternately to obtain improved values of /3 and a. 

To improve the accuracy, start with the approximation obtained graphically, 
and improve it by one or more applications of either of the iterative processes 
discussed in Section 4. 

If b is very small, say of order 10 , %«- <z 3 /4 is one root, and the 
remaining roots can be found from (14) by neglecting 6 , dividing through by y, and 
solving the resulting cubic as indicated in Section 2 - it will be ready-normalized in 
the form (3). 
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